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Abstract. I present a brief overview of the wealth of femtoscopic measurements from the past 
two decades of heavy ion experiments. Essentially every conceivable "knob" at our disposal has 
been turned; the response of two-particle correlations to these variations has revealed much about 
the space-momentum substructure of the hot source created in the collisions. I discuss the present 
status of the femtoscopic program and questions which remain, and point to new efforts which aim 
to resolve them. 



The slowly crawling ants will eat our dreams. 

Andre Biafas, musing on words of Andre Breton as they might apply to femtoscopy. 
Go to the ant, thou sluggard; consider her ways, and be wise. - Proverbs vi.6 



WHEREFORE 

High energy collisions between electrons, hadrons, or nuclei produce highly nontrivial 
systems. Especially in the soft (low-p^, long spatial scale) sector, the inclusive distri- 
butions of the measured multiparticle final states are dominated by phase-space; to first 
order the momentum spectra and particle yields appear thermal, revealing little of the un- 
derlying physics of interest. Detailed information in this sector is obtained only through 
correlations; inclusive spectra tell much less than half the story. 

In particular, multiparticle production is a dynamic process, evolving in space and 
time. For several decades now, small relative momentum two-particle correlations have 
been used to probe the space-time structure of systems at the femtometer scale. Measure- 
ments and constantly-improving techniques variously called "intensity interferometry," 
"HBT," "GGLP," "non-identical correlations," etc, are nowadays discussed under the 
common rubric of femtoscopy [1], as the title of this new workshop series reflects. 

While understanding the space-time features of the system is important to both the 
particle and the heavy ion physicist, in the latter case it is even vital. After all, non- 
trivial geometrical effects dominate the physics of heavy ion collisions. 

From the very broadest perspective, the entire heavy ion program is geared to generate 
and study a qualitative change in the geometric substructure of the hot system. Strongly- 
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coupled |0] or not, the quark-gluon plasma (QGP) is a soft QCD system, in which 
colored degrees of freedom are relevant over large length scales. Of particular interest 
is the existence and nature of a deconfinement phase transition; a significant and sudden 
change in the degrees of freedom should be reflected in space-time aspects of the 
system |0]. Also of generic importance is the (often unasked) question of whether the 
"system" generated is, indeed, a system. Any discussion of "matter" or "bulk" properties 
relies on an affirmative answer. 

More specifically, geometry defines each stage of the system's evolution. In the initial 
state, the entrance-channel geometry (impact parameter b) determines the subsequent 
collective evolution and anisotropic expansion of the system [4]; the resulting "elliptic 
flow" [5] has been the basis of ~ 25% of the publications from the RHIC program. In 
the intermediate state, also, geometry dominates: quantitative understanding of exciting 
parton energy loss (or "jet quenching") measurements 10,0] requires detailed informa- 
tion of the evolving size and anisotropic shape of the system. If coalescence is indeed the 
mechanism of bulk hadronization [ 8 ] , space-momentum correlations in the intermediate 
stage induce clustering effects which must be modeled quantitatively [9]. 

Clearly then, for the soft (bulk) sector in heavy ion collisions, geometrical issues dom- 
inate both the physics of interest and the system with which it is probed. No surprise, 
then, that since the relativistic heavy ion program began roughly two decades ago, fem- 
toscopic studies have played a major role, and a "sub-community" has developed. It was 
not long before the erstwhile "nuclear" physicists contributed physical and technical in- 
sights to a type of measurement initially borrowed from their particle physics colleagues. 
At workshops like this, such dialogue continues unabated. 

Several excellent reviews of femtoscopy in heavy ion physics have very recently 
appeared in the literature Together with physics discussions, the reader 

may find in them precise definitions of the correlation function, "homogeneity lengths," 
"HBT radii," "out-side-long" coordinate system etc. Here, I assume familiarity with 
such concepts, and very briefly review the status of heavy ion femtoscopy at present. 
I emphasize the breadth of systematics which has been explored so far, what (we think) 
it has told us, and what continues to puzzle us. I then identify a few promising directions 
in which the field is moving, pointing for details on these to others' contributions to 
these proceedings. 

WHENCE 

Due to their copious production and ease of detection, most femtoscopic measurements 
have utilized correlations between charged pions. Further, many experiments have fo- 
cused on central (\b\ — 0) collisions, since (1) azimuthal symmetry simplifies the fem- 
toscopic formalism lfl2l fl3ll ; and (2) maximal energy densities and spatial extents are 
generated. The extent of measured femtoscopic systematics 15-20 years ago is repre- 
sented in Figured showing that, in central collisions involving nuclei with mass number 
A, HBT radii scale approximately as A 1 / 3 ifbHHsll . Apparently trivial, these data were 
at the same time comforting, confirming that pion correlations did indeed track with 
geometric scales. 
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FIGURE 1. Pion HBT radius versus the mass FIGURE 2. World dataset of published HBT 
number of colliding nuclei, from Bevalac experi- radii from central Au+Au (Pb+Pb) collisions 
ments ~ 20 years ago. Compilation from flUl . versus collision energy. Compilation from fl2l . 



Since then, femtoscopic data and techniques have evolved tremendously, generating 
an equally tremendous range of systematic femtoscopic studies. The first femtoscopic 
measurement in truly relativistic heavy ion collisions was reported almost twenty years 
ago by the NA35 Collaboration at the CERN SPS Similar measurements have been 
performed at the SPS and the BNL AGS and RHIC accelerators over the collision energy 
range y/sNN ~ 2.3 — 200 GeV. Thus, in each of the two complementary quantities- 
energy and time- we may consider two decades' worth of systematics 11121 . 

The original hope was to find "anomalously" large spatial and/or temporal scales, as 
reflected in the HBT radii, indicating large entropv generation or a long-lived QGP state. 
This expectation was considered rather generic II 1711 . and, guided by quantitative predic- 
tions from hydrodynamical models ll3[ Il8ll . the most commonly-discussed systematic 
was the excitation function (i.e. energy dependence) of pion HBT radii. This is shown 
in FigureEl where no striking features are observed in the HBT radii at any collision en- 
ergy. As I discuss later, this observed contradiction of a seemingly-generic expectation 
may be considered the second "HBT puzzle." 

Clearly, insight into geometrically-driven physics requires more detailed systematic 
studies than the simple excitation function. Indeed, this has always been a generic 
requirement for extracting physics from any observable in heavy ion physics, and has 
required development of heavy ion programs with simultaneous, complementary, large- 
acceptance experiments running at dedicated machines. Especially in the crucial soft 
sector, more is learned by varying independent variables than by long runs at the highest 
possible energy a given machine can deliver. 
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FIGURE 3. Pion HBT radii plotted versus the FIGURE 4. The "effective pion cross-section" 
number of participating nucleons (left panels), and N pio ton ■ Gpx + -Npion ■ Gun an d the "freezeout vol- 
versus the charged particle multiplicity (right pan- U me" ~ Ri ong ■ R 2 side are plotted as a function of 
els). Compilation from (Tj. the collision energy, for central Au+A u ( Pb+Pb) 

collisions. Figure and further details in Kill . 



Inspired by recent "schematic equations" 111 911 . I denote the impressive multi- 
dimensional space explored by femtoscopic experiments as 

Heavy Ion Femtoscopy = R \^/spw',A,B, \b\, 0,)>,my,mi ,m2j (1) 

Global dependences Especially in light of "puzzles," we need to perform a similar study 
as shown in Figure [H checking that femtoscopic radii track with geometric collision 
scales to first order. We may vary the geometric scale of the reaction zone by varying the 
atomic numbers of the colliding nuclei, A and B, and/or by selecting events of varying 
impact parameter, \b\. Of course, fixing only one of these parameters will not define the 
collision scale; instead, a natural quantity would be the number of participating nucleons 
Npart [20]. Pion HBT radii corresponding to different A, B, \b\ and y^vTV are collected 
in Figure The left panels show that these femtoscopic lengths scale similarly to those 
shown in Figure[H replacing A by N part . (Note that results for central collisions, \b\ ~ 0, 
are shown in Figure [H so that N part ~ A.) The HBT radius R ou t, which mixes space 
and time non-trivially, may be expected to violate a pure geometrical scaling; this may 
explain the increased spread in the upper panels of Figure |3] 

To good approximation, at a given a/saw, total multiplicity (a final- state quantity) is 

a function only of N part (an entrance-channel quantity), independent of A, B, or \b\. The 
relationship does, however, depend on collision energy 12211 . As seen in the right panels 
of Figure|3l the final-state multiplicity provides a more common scaling parameter than 
Npart', recent analyses I23LI24LI25|1 show that this scaling persists for different mj values 
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and for lighter colliding systems at RHIC. 

Several observations may be made about this multiplicity scaling. Firstly, it appears 
that knowledge of dN c h/dr\ alone allows "prediction" of the HBT radii (at least Ri on n 
and R s ide)- This suggests that the small increase of these radii with ^Jsnn seen in Figure|2] 
is associated with increased particle production as the collision energy is raised. (Note 
that Np art is approximately constant for the data in Figure El) Secondly, the finite offset 
d in the approximately linear relationship Ri ong ■ R 2 side = c ■ (dN /dr\) +d means that 
freeze-out does not occur at fixed density W24I1 . 

Thirdly, the scaling shown in the figure breaks down dramatically for y/s^N J$ 5 GeV, 
as is obvious from the non-monotonic behaviour seen in Figure |2 As the CERES 
Collaboration has pointed out lElll . this is likely due to the dominance of baryons at 
lower v/^ivw- Indeed, a quantitative connection between the number of protons and pions, 
and a product of HBT radii is possible, by assuming a universal (Ws^- independent) 
mean free path at freezeout Xf. In Figure |H the "freezeout volume" ~ Ri on g ■ Rside 2 
and the "effective pion cross-section" N pio t on ■ c pK +N v i on ■ o KK are seen to coincide by 
scaling the latter by Xf = 1 fm, apparently contradicting the standard assumption that 
freeze-out occurs when the mean free path becomes much larger than the system size. 

HBT radii and the "freeze-out volume" may be connected only the context of a model 
which includes dynamical effects like flow. The analysis of [21] ignores such effects; 
however, its bottom line remains approximately valid, as flow effects on HBT radii are 
expected to be small at low pr II26I1 . 

Kinematic dependences Insight on the dynamical evolution and geometric substructure 
of the emission region is gained by studying the dependence of femtoscopic lengths on 
the next three parameters in "Equation" [TJ 

In non-central collisions, the entrance-channel geometry is naturally anisotropic; the 
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FIGURE 7. Pion HBT radii plotted versus the transverse mass mj for all published measurements of 
central Au+Au (Pb+Pb) collisions over two decades in ^/Itm- Compilation from fT3l . 



hot source geometry approximates the overlap between target and projectile, and is 
characterized by a "long axis" perpendicular to the impact parameter vector b. At RHIC, 
the system expands more rapidly in-plane (|| b) than out (_L b) lEll . If it is, indeed a 
collective system with finite lifetime, then the overall shape should evolve. Pion HBT 

radii have been measured as a function of their azimuthal angle pa i r = Z (k, l^j for 
Au+Au collisions. The measurement at RHIC 112711 is shown in Figure |5j There, it is 
clear that as \b\ — ► 0, the freezeout source becomes larger and rounder. In f act, there 
is a nice "rule of two"- the source expands to twice its original size E H, and 
its anisotropy £ = ((y 2 ) — (x 2 )) / ((y 2 ) + (x 2 )) decreases by the same factor II27I1 . The 
relatively small change in the source shape is at least semi-quantitatively 12911 consistent 
with short timescale estimates [26] based on the longitudinal radius, and at variance with 
expectations from "realistic" simulations IBflll . 

As will become increasingly clear, the only femtoscopic systematic which might 
display non-trivial ^/snn dependence is, in fact, the dependence on <j) p . This is clear 
from the bottom panel of Figure |3 in which the relative paucity of such measurements 
is also clear. It will be especially interesting to see whether the flow and/or timescales at 
the LHC are sufficiently large to produce in-plane freeze-out configurations 

Experiments at a wide range of collision energies have mapped out the rapidity 
dependence of pion HBT radii. Of particular interest here is the so-called Yano-Koonin 
rapidity 7yk IB2L 13311 . which should approximate the rapidity of the fluid element which 
emits a pair of pions at some rapidity Y nK . Figure[6]shows an approximately "universal" 
behaviour YVk ~ Y nK , independent of ■Jsw This i s consistent with (but not proof 
of [12]) emission from a boost-invariant system [33]. 

The most extensively- studied kinematic systematic has been the p^-dependence of 
pion HBT parameters. Figure |7] shows the world dataset of published measurements 
for central Au+Au (Pb+Pb) collisions. The falling dependence of femtoscopic scales 
on transverse velocity is generally believed to arise from collective transverse and 
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longitudinal flow (e.g. $2^). As I mentioned earlier, strong collective flow would be 
an indication that a real bulk system has been formed. 

The longitudinal radius scales approximately as Rj ~ m T 5 , indicating strong longi- 
tudinal flow and again consistent with expectations for emission from a boost-invariant 
system I35LI26I1 . Decreasing transverse radii R () and R s may be due to collective trans- 
verse flow. The simplest flow-dominated models quantitatively interrelate these fem- 
toscopic mj dependences with other observations. An example is shown in Figure 
in which a very simple freeze-out scenario lEr3l - thermal motion superimposed on a 
collectively exploding source- can simultaneously describe a broad range of data mea- 
sured at RHIC. The momentum-space distribution, quantified by the average number 
distribution (top panel of Figure [HJ and the number variation as a function of azimuthal 
angle (middle panel) give an incomplete picture by themselves. Momentum-dependent 
femtoscopic radii (bottom panel) probe the dynamical sub- structure of the collision, 
constraining models more stringently 

Particle-species dependences Within the past several years, high-statistics datasets in 
experiments with good particle identification have allowed the mapping of femtoscopic 
systematics with the final variables in "Equation" [l]- the mass (or species) of the corre- 
lated particles. 

Signals of a system's collectivity at freeze-out should not be limited to the pions. In 
the simplest picture, corresponding to flow-dominated models (e.g. [26]) of Figure [HI 
femtoscopic radii should approximately scale with mj, independent of particle type. An 
impressive common scaling of radii from all measured particles is, indeed, observed at 
all energies explored, as seen in Figure The common scaling is particularly striking 
when one considers the quite different measurement systematics involved in charged 
pion correlations and, say, K® correlations. Even generalized nucleon separation scales, 
probed by relative yields of deuterons and protons (d/p in Figure|9l), follow the system- 
atic, with the exception of one outlier point at the lowest energies. 



TABLE 1. A very incomplete table of published or ongoing femtoscopic studies at RHIC for 
various particle combinations. "Traditional" identical-particle interferometry lies along the lowest 
diagonal line of cells. 
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Correlations between non-identical particles probe not only the sizes, but also the 
relative displacement of the particles' emission zones in space-time W47I1 . Any collective 
freeze-out scenario naturally implies a specific relationship between emission regions 
of the various particle types. In a flow-dominated picture, the emission zones for high- 
niT particles are not only smaller than those for low-mr particles, but are also inevitably 
located further from the center of the collision region 12611 . as suggested by the schematic 
in Figure [10| As discussed in detail in the contribution of A. Kisiel 140I1 . available 
measurements of these displacements at RHIC provide further support of the flow- 
dominated freeze-out scenario. 

Non-identical particle correlations are today a growth industry. Table [T] lists only a 
sampling of recently-published or ongoin g an alyses at RHIC energies. Similar studies 
have been performed at lower energies ll48l [l2ll . The diagonal axis corresponds to 
identical-particle correlations, the "traditional" focus of HBT interferometry. 



WHITHER 

Excluding aficionados attending workshops such as this one, in the minds of most heavy 
ion (or high energy) physicists, the term "femtoscopy" (or, more likely "HBT") brings 
to mind only the single, rather uninspiring systematic plotted in Figure El indeed, some 
may be tempted toward the dismissive view that "the measured radius is always 5 fm." 

As we have just discussed, this is grossly unfair: the systematics are tremendously 
richer, with femtoscopic length scales varying with almost every parameter in "Equa- 
tion" d Furthermore, these strong systematic trends are found consistently by experi- 
ments separated by decades and using quite different measurement and correction tech- 
niques; indeed, especially at RHIC, the data are almost embarrassingly consistent (c.f. 
Figure 13). Yet further, it appears that these systematics may be well understood in the 
commonly-accepted framework of system evolution due to strong flow quantitatively 
consistent with momentum-space observables 12611 . Clearly, there is much more to fern- 




FIGURE 10. Freezeout regions for particles of FIGURE 11. Multiplicity density per participant 
different species (or different transverse masses) pair, measured for Au+Au (Pb+Pb) collisions at 
emitted from a common source. Two-particle AGS, SPS, and RHIC are shown in the lowerpanel, 
correlations measure the (momentum-dependent) taken from 1 22], and naively extrapolated by the 
size, shape, and orientation of the emission re- author to LHC energy, 
gions, as well as the average displacement (Ar) in 
the outward direction. From lTl2ll . 

toscopy than its most notorious Figure. 

On the other hand, the trends shown in Figures |7J and |9] suggest that the 
notorious Figure El quite correctly summarizes the situation after all. At any -Jsnn, the 
systematics of "Equation" [T] are quite rich and may well be reconciled with a reasonable 
physical picture. However, in more ways than expressed by Figure El those systematics 
are essentially independent of ^/snn\ Without resorting to agreement or disagreement 
with particular models, this second 2 femtoscopic puzzle is startling, suggesting that the 
space-time consequences of the physical processes are the same at RHIC as they are near 
the pion production threshold. Often, "universal" behaviour is a key to deeper physical 
insight. Heavy ion femtoscopy, however, might display a bit too much universality. 
Whither... or wither? 



2 In this experimental overview, I have not discussed what has come to be known as "the" HBT puzzle l49ll 
which, simply put, is that otherwise-successful and apparently reasonable models like hydrodynamics do 
not reproduce femtoscopic measurements 1131 . To all but the novice heavy ion physicist, however, the 
initial failure of dynamical models to reproduce diverse observations is hardly puzzling. The experience 
at lower energies is that such initial failure is more the rule than the exception. In light of other, more 
generic puzzles, I call this problem only the "first HBT puzzle" 13711 . 



Given the prominence of nontrivial geometry to the physics of heavy ion collisions 
in general, and the rather generic expectation of significant changes in spacetime 
evolution with ^/snn, understanding this universality remains urgent. What future efforts 
might shed some light? 

Today, one almost reflexively points to the impending heavy ion program at the LHC 
for new observations generating fresh insights. While anything might happen in an un- 
explored energy domain, we may venture a prediction. Instead of a crystal ball, however, 
we use a mirror to gaze over our shoulder at twenty years of systematics in heavy ions. 
The ^/sMv-dependence of the global multiplicity (per participant pair) has been signif- 
icantly extended at RHIC and summarized by the PHOBOS collaboration H22I1 . Boldly 
and probably naively extending this systematic leads to the expectation that multiplici- 
ties at the LHC will be ~ 60% higher than they are at RHIC, as shown in Figure [TT] 

As discussed in the previous Section, femtoscopic length scales- for any mr, y, N par t, 
or particle species- depend primarily on event multiplicity. Taken together, Figures [TTJ 
and|3]suggest that radii 3 in central collisions at RHIC will simply be ~ 17% higher than 
they are at RHIC ((1.6) 1 / 3 = 1.17.). 

Notably, evidence is mounting that perhaps all soft-sector observables are determined 
primarily by total multipilcity, independent of y/spjN- Properly-scaled elliptic flow tlsoh 
and even strangeness enhancement 115 ill appear to show universal multiplicity scalings. 
Whether this is a trivial implication of entropy-driven phasespace dominance in observ- 
ables sensitive to bulk medium is unclear. However, nontrivial new phases of matter 
should have signatures in the long-distance (soft momentum) sector; dependence only 
on multiplicity and not reaction energy would be intriguing. 

So, perhaps the choice of collider facility (LHC versus RHIC) is unimportant, and 
heavy ion femtoscopists should focus on filling in the holes of Table [JJ? Most evidence 
thus far indicates that flow-dominated freezeout scenarios (e.g. 112 611 ) fitted to identical 
K correlations essentially "predict" femtoscopic data using other particle combinations. 
The data is yet scant, however, so ongoing studies Il40ll t o further explore this Table are 
quite important. There are even preliminary reports 113 811 . with exotic particle combina- 
tions, of inconsistencies with these freezeout models. If confirmed, strong theoretical 
focus should come to bear on this result. If, on the other hand, varying the particle com- 
bination repeatedly yields results "predicted" by blast-wave models, continually filling 
in cells of Table [JJrisks becoming a stamp-collecting exercise. 

The most important recent experimental developments in heavy ion femtoscopy were 
presented at these workshops. In addition to those discussed above, I briefly mention 
here some ongoing studies which I find most promising. 

Even if all of the particle combinations in Table [JJ follow simple blast- wave calcula- 
tions, and so reveal no new femtoscopic information, this can actually be turned to good 
use. In particular, one may turn around the traditional approach in which one uses the 
known the two-particle final state interaction (FSI) to extract geometric information, to 
extract the FSI itself HB3l- Finalized results from STAR on p - A correlations IB9ll 



Precise expectations R out or -dependent radii at the LHC are, admittedly, less certain, as the former 
does not scale exactly with multiplicity (cf Figure[3}, and the multiplicity-dependence of the latter has not 
been extensively mapped (cf Figures[2]and[5}- 



have extracted previously inaccessible phase shift information for low-energy baryon- 
antibaryon scattering. While not QGP-related physics, such studies can make a unique 
contribution to low-energy QCD and hadronic physics. 

Ideally suited workshops like these are studies which directly compare for the first 
time, at a fixed energy and using identical detector and analysis techniques, correlation 
data from the heaviest ion collisions to that from p+p collisions [24]. As has been ob- 
served previously in high energy experiments, femtoscopic radii from identical pion cor- 
relations measured in p+p collisions decrease with increasing pj, qualitatively similar to 
the dependence shown in Figure|7J The preliminary STAR data shows, however, that in 
all three HBT radii, the pj dependence is quantitatively identical in p+p and A+A col- 
lisions! Since the heavy ion and high energy communities have traditionally used very 
different physics mechanisms to explain this dependence, this observation potentially 
throws the explanations of both communities into doubt. If this result is confirmed, it 
ranks as the "third HBT puzzle" [37]. 

Unexplained long-range structure in the correlation functions for the lowest- 
multiplicity collisions, however, presently cloud the interpretation of the HBT radii 
for p+p collisions 112411 . Partly in an effort to understand this, a new representation of 
the data in terms of spherical harmonic amplitudes in q- space was developed as an 
experimental diagnostic tool [37]. In fact, a similar harmonic decomposition method 
was earlier already developed by Danielewicz and collaborators J53l l54ll . not as a 
diagnostic, but as a direct link to the detailed geometry (beyond simply length scales) of 
the emitting source. This representation has a natural connection to source imaging Il55ll . 
and, indeed, first applications to PHENIX data have been reported iB^l . 

Harmonic decompositions as an improved representation of the correlation function 
and source imaging as an improved, generalized fit to the data are, in a sense, merely 
technical improvements, but they are quite significant ones. Just as femtoscopic studies 
have explored the systematic landscape of Equation \\\ so should they probe the "mi- 
croscape" of fine details of the measured data. 

The femtoscopy of heavy ion collisions can be an addicting endevour. Systematics 
make sufficient sense that we are convinced that we are probe geometry at the femtome- 
ter scale. Spacetime geometry at that scale is sufficiently important to the physics that 
the measurements must be done well. Such measurements are sufficiently challenging 
that it is enjoyable to do them well and to develop improved techniques. However, for 
now, the overall results are sufficiently puzzling that there is plenty more to do. 
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